Physical implementation of topologically decoherence-protected superconducting qubits 
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We propose a scenario to physically implement a kind of topologically decoherence-protected qubit using 
superconducting devices coupled to a micro- wave cavity mode with unconventional geometric operations. It is 
shown that the two needed interactions for selective devices, which are required for implementing such protected 
qubits, as well as single-qubit gates, can be achieved by using the external magnetic flux. The easy combina- 
tion of individual addressing with the many-device setup proposed in the system presents a distinct merit in 
comparison with the implementation of topologically protected qubits in a trapped-ion system. 
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Physical implementation of quantum computers has at- 
tracted much attention as they are generally believed to be 
capable of solving diverse classes of hard problems. Sys- 
tems suitable for hardware implementation of quantum com- 
puters should possess certain properties, such as relatively 
long coherent time, easy manipulation and good scalability. 
With highly developed fabrication techniques, superconduct- 
ing quantum interference devices (SQUIDs) have shown their 
competence in implementing the qubits for scalable quan- 
tum computation Furthermore, the idea of placing the 
SQUIDs inside a cavity, i.e., the circuit quantum electrody- 
namics, has been illustrated 10. S 0, Btl to have several practi- 
cal advantages including strong coupling strength, immunity 
to noises, and suppression of spontaneous emission. 

Decoherence and systematic errors always occur in real 
quantum systems and therefore stand in the way of physi- 
cal implementation of quantum computers. Generally, larger 
systems are more sensitive to decoherence, which makes the 
scaling of quantum computers a great experimental challenge. 
Therefore, how to suppress the infamous decoherence effects 
is a main task for scalable quantum computation. A promis- 
ing strategy for quantum computation in a fault-tolerant way 
is based on the topological idea where gate operations 
depend only on global features of the control process, and 
are therefore largely insensitive to local noises. The Kitaev 
model 1 6] consists of a class of stabilizer operators associated 
with lattice on the torus, which can be put together to make 
up a Hamiltonian with local interactions. The four-fold de- 
generate ground state of the Hamiltonian coincides with the 
protected space of the stabilizer quantum code. Since all the 
excited states are separated from the ground states by an en- 
ergy gap, the ground states are persistent to local perturba- 
tions. However, it is extremely difficult to directly implement 
this novel idea mainly because four-body interactions are no- 
toriously hard to generate in experiments. 

Very recently, Milman et al. [TJ proposed a highly symmet- 
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rical Hamiltonian 

M ( M \ ^ M / M \ ^ 

which involves only two-body interactions and can be under- 
stood as spins in a 2-dimensional (2D) lattice. Each spin was 
labeled by its position in the lattice with a^'^ denoting Pauli 
matrices of the spin situated at the intersection of the ith row 
and jth column. It has a two-fold degenerate ground state 
which may be considered as the two states in a qubit. These 
states are generated by non-local operators and thus are nat- 
urally protected against local sources of decoherence. How- 
ever, comparing with a size-independent energy gap in the Ki- 
taev Hamiltonian, the energy gap here is slightly dependent on 
the number of the involved spins. A nice feature of the model 
lies in that it is possible to directly implement the Hamilto- 
nian ([T]i in a trapped-ion system, as it eliminates a main ex- 
perimental constrain-the four-body interactions in the Kitaev 
model. In addition, the energy gap of Hamiltonian ([T]) remains 
large as the system size increases, and thus it provides more 
efficient protection than the nearest neighbor model 1 8] where 
energy gap decreases rapidly. 

In this paper, we propose an intriguing scenario to imple- 
ment the model Hamiltonian ([T]i with superconducting devices 
(serve as spins in ||^01) coupled to a micro-wave cavity mode 
with unconventional geometric operations. In trapped-ion sys- 
tems, the individual addressing is difficult for large arrays 
{N > 3) since the distance between the ions in the center of 
the linear array gets smaller ||^[l3l- This constrain makes the 
linear configuration proposed in Ref. |01 difficult to be experi- 
mentally implemented. Comparing with the above difficulty, a 
distinct feature of the present implementation is that the com- 
bination of individual addressing with a many-device setup is 
feasible. We will show that, by using the external magnetic 
flux as the effective switch tool, the interactions between se- 
lective devices can be introduced. Furthermore, single-qubit 
gates, which are required in the global state initialization, can 
also be achieved by the same technology. During the imple- 
mentation of the operations, the only parameter needs to be 
tunable is the frequency of the external magnetic flux. This 
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FIG. 1: (Color online) Schematic illustration of the superconduct- 
ing device as the effective spin, the red dash line denote the integral 
path of the magnetic flux, (a) A single SQUID as an effective spin, 
(b) Device made of two SQUIDs with a common superconducting 
charge box. This more flexible design will introduce more control 
variables of the effective spin. 



is a more efficient way of controlling the system dynamics 
as it is much easier for experimental realization. In achiev- 
ing the implementation of the needed interactions and given 
the fact that selective coupling can be effectively controlled, 
it is possible to implement topologically protected qubits in 
this architecture. Moreover, the operations to achieve the 
above qubit states are based on unconventional geometric op- 
erations, which have been illustrated to have a high fidelity 

A SQUID as a superconducting charge qubit that usually 
addressed in the literature is shown in Fig. ([T^) How- 
ever, it would be clear that this simple configuration can not 
be used to acheive the coupling described by Eq. ([T]). To 
realize the desired Hamiltonian, we design a more flexible de- 
vice, as shown in Fig. dTJi), to serve as the effective spin. 
It consists of two SQUIDs with a common superconducting 
charge box, which has n excess Cooper-pair charges. Each 
SQUID is formed by two small identical Josephson junctions 
(JJs) with the capacitance Cj and Josephson coupling energy 
Ej, pierced by an external magnetic flux. A control gate volt- 
age Vg is connected to the system via a gate capacitor Cg. The 
Hamiltonian of the system reads [IJ 



H = Ecin - nf 



E 



E 
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cos Lp' 
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where n is the number operator of (excess) Cooper-pair 
charges on the box, Ec — ^e^ /{Cg + 4(7j) is the charging 
energy, n = CgVg/ 2 is the induced charge controlled by the 
gate voltage Vg, and tp^ (I = \,2, 3, 4) is the gauge-invariant 
phase difference between the two sides of the Zth JJ denoted 
as Ji in Fig. ([TJ)). The phase differences ip^ are determined 
from the flux quantization for three independent loops, that 
is, ip'^ ~ ip^ ^ - ^3 ^ 20, and ^3 - ^ 20^. 

Since we here focus on the charge regime, a convenient basis 
is formed by the charge states, parameterized by the number 



of Cooper pairs n on the box with its conjugate Lp = J^if'' 
by the relation [ip, n] = i. At temperatures much lower than 
the charging energy and restricting the gate charge to the range 
of n G [0, 1], only a pair of adjacent charge states {|0), |1)} 
on the island are relevant. Under the condition (p^ — (p'^ ~ 0, 
the Hamiltonian ^ is then reduced to 



-Ececr'' - 2E^a'' 



(3) 



where Ece = 2Ec{l ~ In), E^j, = Ej cos {n^/tpo) = 
Ej cos (j) with (po — hc/e being the normal flux quantum and 

(j) = TT<^/(j)Q. 

When the device is placed in a cavity, the gauge-invariant 
phase difference becomes (p,,, = ipm ^ ^ Ji ■ dim, 

where A„i is the vector potential in the same gauge of ipm- 
A„i may be divided into two parts A'^^ + A^^, where the 
first and second terms arise from the electromagnetic field 
of the cavity normal modes and the external magnetic flux, 
respectively. For simplicity, we here assume that the cavity 
has only a single mode to play a role. Therefore, we have 



2tt r 



2-n- r 
00 -llj 



dl„ 



2g\ 



where 



2g is the coupling constant between the junctions and the cav- 
ity, a (flt) is the annihilation (creation) operator for the cavity 
mode, with loc being its frequency, the closed path [the red 
dash line in Fig ([T]l] integral of the A"^ gives rise to the mag- 
netic flux Then the Hamiltonian ^ becomes 



He = -Ece(Jz - 2Ej cos y> + g (a + a^)] 



(4) 



We now show that the required single-qubit gates can be 
achieved by the dc magnetic flux. Generally speaking, since 
the coupling constant g is very small comparing with 0, we 
may expand the Hamiltonian ^ up to the first order of g. 
In the rotating frame with respect to Hdc ~ —EceO'z — 
2E^c7x + f^c {o)a + i) , the strength of the cavity-SQUID 
interaction term in Eq. Q is proportion to l/(wc ~ ^q) with 
ijjq = 2Ec{l — 2n)/h, which is an experimentally controllable 
parameter via the gate voltage Vg. In the present scheme, we 
choose LUq = 0, which corresponds to the degeneracy point. 
Given the facts that g is relatively small and the frequency dif- 
ference LUc — is chosen to satisfy the large-detuning limit, 
the cavity mediated interaction can thus be safely neglected. 
In other words, the device and the cavity evolve independently 
in this case. The external flux is merely used to address sep- 
arately the single-qubit rotations, while the evolution of the 
SQUID is governed by Hamiltonian (O. 

The required interactions between the selected devices may 
be induced by using the ac magnetic flux with devices in Fig. 
([TJ)). If N devices are located within a single-mode cavity 
as shown in Fig. (|2^), to a good approximation, the whole 
system can be considered as A'^ two-level systems coupled to 
a quantum harmonic oscillator fs*]. The Hamiltonian of the 
such devices, placed in a single mode cavity, reads ||5il 



N 



Eceirij 



^2 - Ej^ COS (p\ 



I 



(5) 



where we have assumed Ece.j = Ec. 



and E'j .j = Ej. Here 



and (f> are the dc magnetic flux, and we have neglected the 
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FIG. 2: (Color online) (a) Schematic superconducting devices (de- 
noted as rectangles) array constructed along the cavity direction x 
with A'^ = 9. The chosen coordinate is shown in the left-lower panel. 
Each device is placed parallel to the xoz plane and at the antinodes 
of the single-mode standing-wave cavity. The magnetic component 
of the cavity mode is along the y direction, which is perpendicular 
to the device loop plane, so that it is the only contributed compo- 
nent. Devices are placed at the antinodes of the cavity mode, so that 
the device-cavity interaction constant {2g) of different devices can 
be treated as the same constant, (b) An A'^ = 9 linear array corre- 
sponds to a 2D 3x3 array, where the filled dots with the number 
represents the devices. For example, if one wants to generate the in- 
teraction in row 1, then devices 1, 2 and 3 are selected to achieve the 
coupling via the virtue excitation of the cavity mode, while the 
Jy coupling among devices 2, 5 and 8 is implemented in order to 
achieve the interaction in column 2. 



cavity-mediated interaction parts in these two SQUID loops 
considering that they can be made to be very small comparing 
with the inter-SQUID loop, or we can just put the magnetic 
shield to exclude the influence of the cavity mode magnetic 
flux in these two loops. The phase relation in the inter-SQUID 
loop is now modified as (p'j — (pj = 2(j)j + 2gj{a + a^) = 
2u!t + 2g{a + a^). Assuming all the devices work in their de- 
generacy points, under Lamb-Dicke (LD) limit and rotating- 
wave approximation as well as in the interaction picture with 
respect to iJo — ?iijJc(a^a+|), the cavity mediated interaction 
can be described by the Hamiltonian 



N 



i5t\ 



H.C., 



(6) 



where gEj/2h ^P^6 = lUc-lu-^uj, a"^^ i((T^ iicr?'), 

and = kn, then Eq. (|6]l 



reduces to 



(7) 



Hamiltonian is given by Hill 

where Xx = 4/3^/(5. If 0_ = 7r/2 and 
Eq. ^ reduces to 



(8) 

k-K + tt/2, then 
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where X 



J2j=i fj'^'^- The corresponding effective 



2ihp [a^e 

with the effective Hamiltonian being given by 

Hy = -hXyJy, 

where Xy — Xx- Actually the configuration described in Fig. 
([l^) is exactly equivalent to the device in Fig. ([TJ)) in the 
specific case of cj)^ — cj)^ = 0. So it is now seen that the 
required Hamiltonian ([T]| is unable to be achieved with the 
one SQUID in Fig. ([T^) since only the interaction Hx could 
be realized. 

It is clear that the specific choice of the dc magnetic flux can 
lead to the designated interaction type of the selected device. 
The cavity-device coupling and decoupling can be controlled 
by selecting the external magnetic flux to be dc or ac, not sim- 
ply by changing the parameters of the device or the cavity. 
Since the external magnetic flux can be effectively controlled, 
the cavity-device interaction can be implement in the selected 
devices. In addition, all the devices are always stay in their 
degeneracy points in this implementation, where they possess 
long coherence time and minimal charge noises. Interestingly, 
the interactions in Eqs. (O and ( fTOl i. which correspond to the 
row and column interactions in Eq. ([T]), are insensitive to the 
thermal state of the cavity mode. It is notable that the evo- 
lutions governed by the Hamiltonians ([8]) and ( fTOl i may also 
be considered as unconventional geometric operations ITTIl . 
which are robust against random operation errors ifisll . and 
thus have been essent ially used in quantum information pro- 
cessing i, dm in III [H. 

With the two interactions in Eqs. ([8]) and ( fTOb . and given 
the fact that they can be implemented on selected devices, it is 
straight forward to implement the highly symmetrical Hamil- 
tonian ([T]i. The interaction exists only between the devices 
selected at one time, thus the spatial pattern of the physical ar- 
ray is not necessarily related to the 2D configuration directly. 
A possible configuration rendering the implementation of the 
Hamiltonian easier in existing systems consists of using a lin- 
ear array of devices [7]. To implement a topologically pro- 
tected logical qubit consists of AI x M physical devices in 2D 
lattice, we can use a ID linear array with N = devices, 
where the ith row and jth column in the 2D lattice correspond 
to the devices (i - 1) x M + 1, (i - 1) x M +2,...,ix AI and 
j, j + M, j + (M — 1) X M in the ID array, respectively. 
A small 3x3 2D array is illustrated in Fig ^p) as a example. 
After scaling the interactions of each row and column, the re- 
sulting Hamiltonian is the sum of the terms in Eq. ([T]i. The 
coupling between rows and columns may be avoided |7] by 
applying operations that alternate between rows and columns 
provided that each operation time r satisfies TXx;y <C 1. 

Topologically protected qubits usually consist of many 
physical qubits and thus seem to be resource-consuming strat- 
egy. In fact, a small finite array, e.g., M = 3 for Hamilto- 
nian and M — 5 for nearest neighbor model ISt], are 
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likely sufficient to provide good protection from the noise by 
suppressing its effect to be many orders of magnitude lower, 
which could meet most practical demands. 

One experimentally feasible procedure to generate a topo- 
logically protected qubit state is as follows, (i) Applying a 
large effective field along the x direction (putting n = 1/2 
and (f) — Q or (f) ~ it) when devices are decoupled from the 
cavity, we get Hi ~ ±Ej J2 '^fj ■ Then the initial state is log- 
ical |0) = |0)y or logical |1) = H depending on the 
overall sign in the Hamiltonian Hi. (ii) Adiabatically switch- 
ing off Hi and then switching on the Hamiltonian ([T]). After 
performing these two steps, the two initial unprotected logi- 
cal states |0) and |1) evolve into one of the protected ground 
states of the model Hamiltonian ([T]i. 

We now briefly address the experimental feasibility in im- 
plementing our scheme. Suppose that the quality factor of 
the superconducting cavity is Q = 1 x 10^. For the cav- 
ity mode with lUc/^tt — 50 GHz [17], the cavity decay time 
is Tc « 3.2 /IS. With a moderate vacuum Rabi frequency 

- 10 MHz and the lifetime of the SQUID 7 = 2 //s, the 
strong-coupling limit can be readily reached [4]. Meanwhile, 
2g = 10"^ 10] to ensure the LD limit. With Ej = 40 //eV 
yD, /3/27r « 48 MHz. To maintain the large-detuning condi- 
tion, we can choose S ~ 10/3, which in turn satisfies readily 
the requirement of S -^i oj ^ ojc- A typical operation time 
for large-detuning interaction is t w 10 ns |5], which ensures 
that thousands of operations are possible [4] since the coher- 
ence time of the SQUID and cavity mode are both in the or- 
der of 1.1s. When all devices are located at the antinodes of 
the microwave, the distance of the two neighbor devices is 
a half of the microwave wavelength. For the chosen cavity 



mode (Ac = 6 mm), 10 devices may be constructed along 
the cavity direction iflTll . TypicaUy, a SQUID loop (1x1 
/im) consists of two small identical JJs of the size 0.1 x 0.1 
/im, and the inter-SQUID loop can be fabricated in the size 
of 4 X 4 iim. Therefore, the distance of the two devices is 
about 10'^ times of the device size, so that their mutual in- 
duction may be safely neglected. Moreover, it is notable that 
other kinds of 4-junction qubits have already fabricated and 
used in the experiments, e.g., in jlStl . Thus, the current design 
of qubit impose no extra requirements on the current tech- 
nology of the qubit fabrication process. Finally, we roughly 
estimate the effect of the nonuniformity in device parameters, 
Scejand Ej, to the operations in our implementation. As dif- 
ferent devices can be addressed and controlled individually, it 
is easy to tune each device in its degeneracy point via its gate 
voltage separately, such that the nonuniformity of Ece among 
different devices is not of significant importance. The nonuni- 
formity of the Josephson couplings with a deviation of eEj 
will cause a minor change of the operation infidelity in the or- 
der of exp{—A/£Xx), which is about 1% and thus can be ne- 
glected for the reasonable experimental parameters e ~ 20% 
and S ~ 10/3, where A ~ is the energy gap of topological 
phase. 
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